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Abstract 

The k±_ dependent gluon distribution is calculated numerically accounting for the resummation of 
small x effects due to the Lipatov equation. It is represented by a convolution of a gluon density in 
the collinear limit and a universal function Q(x,k 2 ,fi) for which an analytic expression is derived. 



1. Introduction 

In the small x range new dynamical effects are expected 
to determine the behaviour of structure functions. The 
evolution of parton densities is effected by terms due 
to non strong k±_ ordering and eventually by screening 
contributions. A description of contributions of this 
type requires to generalize the factorization of the 
hadronic matrix elements into coefficient functions and 
parton densities in which the k±_ depencence is not 
integrated out [1]. 

This factorization covers the case of collinear 
factorization in the limit that the k 2 dependence of 
the coefficient function is neglected. The k±_ dependent 
gluon density accounts for the resummation of small 
x effects. In the present paper we will consider those 
due to the Lipatov equation only. Since this equation 
behaves infrared finite no other singularities will emerge 
than in the case of mass factorization. The collinear 
singularities are delt with in the same way in the case 
of k±_ factorization. 

In the present paper the k±_ dependent gluon 
distribution is calculated for the case of the scheme 
[2, 3]. It can be represented as the convolution of 
the gluon density in the collinear limit g(x,/i) and a 
function Q(x,k 2 ,fi) for which an analytic expression 
will be derived. A numerical comparison of the k±_ 
dependent gluon densities accounting for the solution 
of the Lipatov equation and in the double logarithmic 
approximation (DLA) is given. 

* Talk given in session on Proton Structure at the Workshop on 
Deep Inelastic scattering and QCD, Paris, April 1995 



2. k±_ Factorization and the k±_ dependent 
gluon distribution 

The factorization relation for an observable Oi(x,/i) 
reads 

Oi(x,fx) = J dk 2 a 0i {x,k 2 ,ii)®§{x,k 2 ,ii) (1) 

where aOi(x,k 2 , fi) and Q(x,k 2 , fi) denote the k 2 
dependent coefficient function and parton density)), 
respectively. Eq. (1) can be rewritten as [2, 3] 

Oi(x,fi) = o-%.(x,fi)®G(x,fi) 

+ J^dk 2 [a Oi {x,k 2 ,fx)-a- Oi {x,fx)]^{x,k 2 ,fx) { > 

with aQ,(x,fj,) = lim.i.2^Q aOi(x, k 2 , fi). The first 
addend in (2) describes the conventional contribution 
due to collinear factorization. The second term contains 
the new contributions. Note that Q(x,k 2 , fi) starts with 
terms oc a s . It has therefore not the interpretation of a 
probability density and may even become negative. 

As shown in [2] the k±_ dependent gluon distribution 
associated to eq. (2) reads in moment space 

1 / , k 2 V c{j '"' ) _ 
<S>{j,k 2 ,fi)= 7c {j,a s )—l — \ g{j,fx) (3) 

where /i denotes a factorization scale, a s = 7\T c a s (/i)/7r, 
and g(x, /i) is the gluon density. Eq. (3) accounts for the 
small x behaviour due to the Lipatov equation. Here 

[j We will consider the gluon density in the present paper only. 
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Fig. 1 Real and imaginary part of y c vs i?e(p) and 7m(p). 
In x space the k±_ dependent distribution is given by 

ffln several recent approaches [4] to describe kj_ dependent 
gluon distributions phenomenological Ansatze were used based on 
solutions of inhomogeneous Lipatov equations. Note that these 
descriptions are not related to eq. (2) and [2, 3]. 



<P(x,k ,fi) = y(x,k ,fi)®g(x,fi), 
correspondingly, with 



Jo 



dk 2 $(x,k 2 ,fi) = 8(1 - x). 



(5) 



(6) 



The function Q(x,k 2 ,fi) is universal and can be 
calculated numerically by a contour integral in the 
complex plane over the first factor in eq. (3). Since the 
solution of eq. (4) is multivalued the Mellin inversion 
to x space requires to select the branch in which 
for asymptotic values of j e C j c approaches the 
perturbative result j c (j, «s) ~ &s/(j — 1) for small 
values of a s . 

We solved eq. (4) under this condition numerically 
using an adaptive Newton algorithm. The solution is 
characterized by the well-known branch point at p = 
(j — l)/a s = 4 In 2 and two further conjugate branch 
points (cf. also [5]). As shown in figure 1 the 'ridge' 
in the real part due to the singularity at p = 4 In 2 
persists until Rep ~ —1.3 and turns into a flat form 
with Re(j c ) ~ —0.5 for |/mp| < 1.5. At the same time 
the imaginary part of j c becomes continuous again. 

3. An analytical solution for Q(x,k 2 ,p,) 

Because the integration contour has to be situated 
outside the range of the singularities of j c one may 
expand j c (j, «s) into a Laurent series over p 



7c (j, a,) 
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The coefficients gi are given in [6] up to / = 20 in 
analytical form extending an earlier result [7] . For small 
values of |/mp| (|Jmp| < 2) the truncated Laurent 
series leads to an oscillatory behaviour and eq. (7) does 
no longer serve to be an appropriate description of j c 
(cf. [6]). 

Using (7) a corresponding expansion may be 
performed for 

k 2 G{j,k 2 ,fi) = j c (j, a s )exp[7 c (j, a,)L] (8) 

with L = ln(fc 2 //i 2 ). For the single terms of the 
Laurent series in p the Mellin transform can be carried 
out analytically. Here it is important to expand the 
exponential in eq. (8) in such a way that the lowest 
order term in a s of j c is kept in exponential form. One 
obtains 

k 2 G{x,k 2 ,p) = ^-I (2i/a. log(l/aj)L) 



x !„_! (2 Va.log(l/aj)£) ,£>(). (9) 

The coefficients d v (L) are given in ref. [6]. Up to v = 20 
they contain at most terms oc L 4 . The first term in 
eq. (9) denotes the Green's function in DLA. 
For L — > (9) takes the form 



fc 2 c?(M 2 ,M) = — £• 



9! 
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(10) 



and for L < (i.e. fc 2 < /i 2 ) one has 
k 2 g{x,k 2 ,fx) = ^-J (2^a s log(l/aj)|£|) 

x J„_! (2V«.log(l/a!)|£|) • (11) 

Thus for fc 2 — > damped, oscillating modes are obtained 
which vanish faster than 1/|_L| -1 / 4 . 



result for k 2 ~ /i 2 at a; ~ 10" 3 by 10 to 15% while for 
k 2 — > smaller values are obtained. At larger values of 
x the complete solution approaches the DLA result. For 
k 2 — > $(a:, k 2 , fi) vanishes. Since the DLA result is pro- 
portional to J (2 v / a s log(l/a;) log \k 2 /fi 2 \) for k 2 — > 
a damped oscillatory behaviour is obtained in this ap- 
proximation. The complete solution, on the other hand, 
behaves monotonous in the whole kinematical range. 

The calculation of G(x, k 2 , fi) in eq. (5) by numerical 
Mellin inversion using the numerical solution of eq. (4) 
is rather time consuming compared to the convolution 
of the analytical solution (sect. 3) with the input 
distribution g(x,/i). We compared both methods and 
found that the representation given in the previous 
section leads to a relative error of less than 0.002 using 
an expansion up to 0(a 20 ). 

Since shape and size of the complete solution and 
in DLA are rather similar very precise measurements 
are required to establish the non-DLA contributions at 
small x. 

5. Conclusions 



4. Numerical Results 

The k]_ dependent gluon distribution ^(x,K 2 ,fi) 
(scaled by k 2 ) is shown in figure 2 as a function of x 
and k 2 for fi 2 = 20GeV 2 . 
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Fig. 2 The k± dependent gluon distribution as a function 
of k 2 and x. Full lines: complete solution eq. (5); dashed 
lines: solution in DLA. For the input distribution g{x,fi) 



We have calculated the k±_ dependent gluon density 
numerically in leading order using the BFKL equation. 
A consistent treatment of observables is possible in 
the scheme [2, 3]. The Green's function Q(x,k 2 ,fi) 
was found both numerically and by a perturbative 
analytic expression expanding the complete solution up 
to 0(a 20 ). Both representations agree better than 0.002 
after convoluting with the respective input distributions 
g{x,fi). 

The effect of the non-DLA terms in Q(x,k 2 , fi) is of 
O(10. ..15%). To reveal these contributions requires very 
accurate measurements in the small x range in the case 
of every observable being sensitive to the gluon density. 
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